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Abstract 

The temporal evolution of the combustion process estab- 
lished during projectile transition from the launch tube into 
the ram accelerator section containing an explosive hydro- 
gen-oxygen-argon gas mixture is studied. The Navier-Stokes 
equations for chemically reacting flow are solved in a fully 
coupled manner, using an implicit, time accurate algorithm. 
The solution procedure is based on a spatially second order 
total variation diminishing scheme and a temporally second 
order, variable-step, backward differentiation formula 
method. The hydrogen-oxygen chemistry is modeled with a 
9-species, 19-step mechanism. The accuracy of the solution 
method is first demonstrated by several benchmark calcula- 
tions. Numerical simulations of two ram accelerator configu- 
rations are then presented. In particular, the temporal 
developments of shock-induced combustion and thrust 
forces are followed. Positive thrust is established in both 
cases; however, in one of the ram accelerator configurations 
studied, combustion in the boundary layer enhances its sepa- 
ration, ultimately resulting in unstart. 

Introduction 

The ram accelerator is a chemical propulsion method for 
accelerating projectiles to very high speeds. In this device, 
developed at the University of Washington [1-3], a projectile 
is accelerated inside a tube filled with an explosive gaseous 
mixture, as shown schematically in Fig. 1. The ram accelera- 
tor consists of a light gas gun (using helium as the driver 
gas), ram accelerator section, final dump tank and projectile 
decelerator. The light gas gun provides the initial accelera- 
tion to the projectile, which travels through the evacuated 
launch tube and enters the ram accelerator section. The large 
evacuated tank serves as a dump for the helium driver gas. 
The ram accelerator section can be divided into several sec- 
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tions containing different explosive gas mixtures. The vari- 
ous sections are separated from one another by thin Mylar 
diaphragms, which are ruptured by the projectile. 

Several ram accelerator operation modes, sp anning the 
velocity range 1-12 km/s, have been proposed. In the high- 
speed mode of operation considered in this paper, ignition is 
achieved by means of shock heating. The combustion reac- 
tion couples with the shock, forming either a shock-induced 
combustion wave or a detonation wave, depending on the 
mixture composition, pressure and tube diameter. 

Computational studies of this concept have in the past 
been either focused on steady-state solutions [4-9] or based 
on fully transient but inviscid approaches, utilizing simple 
one-step combustion models [10-12], The first time-accu- 
rate, viscous simulations of the ram accelerator were 
reported by Nusca [13], who used a global methane mecha- 
nism. 

A transient flow calculation clearly provides a more 
accurate description of the combustion process in a ram 
accelerator than does a steady-state solution. Also, it has 
been shown that simple global reaction mechanisms often 
cannot accurately describe shock-induced combustion, 
because of their inability to predict induction times correctly 
[5]. Furthermore, recent calculations [4,5] have demon- 
strated that viscous effects are of p rimar y importance, not 
only in ram accelerators but in any hypersonic propulsion 
concept. These observations illustrate the need for time- 
accurate, viscous (tuibulent) computations using detailed 
combustion mechanisms. 

Such calculations could not hitherto be accomplished, 
due in large part to the lack of an efficient numerical algo- 
rithm. Most previous time-accurate simulations of multi-spe- 
cies reacting flows used explicit or point implicit methods, in 
which only the chemical source term was treated implicitly. 
The reason is partly because the governing equations 
become stiff, thereby complicating the chemical source 
terms for the commonly used implicit methods. Since 
explicit or point implicit methods are constrained by the 
CFL condition, they are very inefficient for solving viscous, 
reacting flows. 

Recently, we developed a computational fluid dynamics 
(CFD) code that addressed the need for efficient time-accu- 
rate simulations of chemically reacting, viscous flows [14], It 



is based on a spatially second order total variation diminish- 
ing (TVD) scheme and a temporally second order, implicit, 
variable-step, backward differentiation formula (BDF) 
method. The inversion of large matrices is avoided by parti- 
tioning the system into reacting and nonreacting parts; a 
fully coupled interaction is, nonetheless, maintained. 

In the present paper we use this code to study the tempo- 
ral evolution of the shock-induced combustion process in a 
ram accelerator. In particular, we investigate the transition 
from the launch tube into the ram accelerator section con- 
taining an explosive hydrogen-oxygen-argon gas mixture. 
The hydrogen-oxygen chemistry is modeled with a 9-spe- 
cies, 19-reaction mechanism [14]. 

Two simplifying assumptions are made in this work, 
which represents a pre liminar y numerical study of flow 
establishment in a ram accelerator. First, the bursting of the 
diaphragm occurs instantaneously and ideally, at the moment 
of projectile arrival. Second, in order to avoid computations 
involving extremely large pressure gradients between the 
evacuated launch tube and the first ram accelerator section, a 
buffer section containing an inert gas is introduced between 
the two segments, as shown in Fig. 2. 

Numerical Formulation 

Governing Equations 

The conservation form of the nonequilibrium Navier- 
Stokes equations describing two-dimensional or axisymmet- 

ric chemically reacting flow involving n species can be 
written in general curvilinear coordinates (^, q ) as follows: 


dt as 


d(G-G v ) 

+ 5n 


+j(s-s v ) 


= w 


( 1 ) 


terms in Eq. (1) and additional state and constitutive equa- 
tions needed for system closure are given by Yungster [4]. 

Numerical Method 


The numerical method used for solving Eq. (1) is 
described in detail in Ref. [14]. Here we present only a brief 
s ummar y of the algorithm. For simplicity, only the two- 
dimensional Euler equations are considered in this descrip- 
tion, although extension to viscous flows is straightforward 
[14]. The equation set is discretized using a temporally sec- 
ond order, variable-step, BDF method. The difference equa- 
tions are solved step by step; that is, approximate solutions 


{Q. .} are generated at the discrete time points t 
J> K 

(n — 1,2,...). Thus, starting with the known initial condi- 
tions {<2° p at t = ?q the numerical method advances the 
solution at each grid point in time, until the desired end state 


is reached. At each integration step, the time step At” is 
selected automatically, by using the procedure described by 
Yungster and Radhakrishnan [14]. 


Ass uming that approximate solutions have been pro- 
duced at the times f” J (j = 0, 1, ...), the formula for 
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advancing the solution to the current time t (i.e., for 
solving the governing equations over the time interval 
T n n + 1"1 , 

1 1 ,t J ) can be written as: 


(3) 

" F j - 1/2, it + G j, k + 1/2 ” G j, it - 1/2 " W j, it ] 


where the parameter j is zero for two-dimensional flow and 
one for axisymmetric flow, and Q is the vector of dependent 
variables: 


where, for the current step. 


. „n+ 1 

AQ ik = Qj,k 


Q 


■j.k 


(4) 
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Q - J [Pj, P 2 P w > P“- P v - e ] (2) 

The dependent variables are the mass density of the z'th 
species p^, the velocity components u and v, and the total 

energy per unit volume e. J is the grid Jacobian, and F and 
G are the inviscid flux vectors in the E, and q directions, 
respectively. Similarly, and G^ are the viscous fluxes. 

The terms S and are axisymmetric source terms, and W 
is the chemical source term. A detailed description of the 


is the incremental solution vector, y and P are the variable- 
step BDF method coefficients [14] and At 1 (= + 1 - t”) 

is the time step. The terms F and G are the numerical 
fluxes in the ^ and q directions. They are computed using 
Yee’s second order TVD scheme [15]. 

Equation (3) is linearized in a conservative manner and 
solved iteratively, using a lower-upper relaxation procedure 
consisting of successive Gauss-Seidel (LU-SGS) sweeps. At 
each time step, successively improved approximate solutions 
to Eq. (3) are generated, until an appropriate convergence 
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criterion is satisfied [14]. For the method to remain tempo- 
rally second order, the flux Jacobians (= dF/dQ and 
dG/ 3(2) must be exact; otherwise, the accuracy reduces to 
first order. The inversion of large matrices is avoided by par- 
titioning the system into reacting and nonreacting parts; 
however, a fully coupled interaction is preserved. The bene- 
fit of the partitioning is that the computational cost of the lin- 
ear algebra associated with matrix inversion is the same as 
that for the commonly used point implicit methods. The rea- 
son is that the matrices arising in the two approaches are of 
the same size. Another important advantage of our method is 
that it remains stable for large values of the CFL numb er, 
thereby enabling the use of large time steps. 

Results 

The goal of this work was to study numerically reacting 
flow establishment during projectile entrance into the ram 
accelerator section. However, before attempting to solve this 
problem, the accuracy of the method was assessed by solv- 
ing various time-dependent flows for which experimental 
data or results of previous numerical simulations were avail- 
able. Two such “benchmark” cases, involving both reacting 
and nonreacting flows, are presented below. 

Benchmark test cases 

The first case was a simulation of Lehr’s [16] ballistic 
range experiments, which consisted of spherical nosed pro- 
jectiles of diameter 15-mm being fired into a premixed, sto- 
ichiometric hydrogen-air mixture. Figure 3a shows the 
shadowgraph image obtained by Lehr [16] for a Mach num- 
ber M = 4.79. The corresponding computational result 
obtained with a 220x220 uniform grid and a 9-species, 19- 
step reaction mechanism for the hydrogen-oxygen chemistry 
[14] is shown in Fig. 3b in the form of density contours. 
Under the conditions of the test, the reacting flow was unsta- 
ble, resulting in a highly regular, periodic flow structure. An 
experimental frequency of oscillation of 720 kHz was 
reported [16]. The computed frequency varied from 701 to 
716 kHz. Computations for other flow conditions also exhib- 
ited excellent agreement with experimental data [14], 

The second test case, taken from the work of Young and 
Yee [17], simulated shock wave diffraction from a 40° 
wedge in air, as shown in Fig. 4a. The top of the wedge was 
rounded, with a radius of curvature of 0.17 times the base 
width. Figure 4 presents the temporal evolution of the air 
flow during the nonreacting diffraction process for an inci- 
dent shock Mach number of 2.0. A 313x140 uniform grid 
was used, and inviscid flow was assumed, after Young and 
Yee [17], Figure 4b shows density contours during the for- 
mation of the triple point, with the Mach stem and a contact 
discontinuity emanating from it. As the shock wave moves 
over the rounded top, the Mach stem evolves into a curved 
shock (Figs. 4c and 4d), which travels slightly faster than the 


incident shock. Figure 4e shows the reflection of the shock 
from the flat rear portion of the body and the formation of a 
new Mach stem. Figure 4f shows the computational results 
obtained by Young and Yee [17] at approximately the same 
instant as that of Fig. 4e. Their calculations were based on 
the MacCormack symmetric TVD method [18]. Hie qualita- 
tive agreement between our calculations and theirs is excel- 
lent The solutions of Young and Yee [17] at previous times 
(not shown here) were also in very good qualitative agree- 
ment with our calculations. 

Ram accelerator reacting flow establishment 

Two projectile configurations are presented. In the first 
case the projectile tail was truncated to resemble closely the 
projectiles used in the University of Washington experiments 
[1-3]. In the second case, the tail ended at a sharp point, and 
the projectile shape was modified with the aim of maximiz- 
ing thrust Both cases modeled the transition from a pure 
oxygen buffer section at o = 1 atm, T = 300 K into a ram 
accelerator section con tainin g an explosive mixture of 
hydrogen, oxygen and argon at the same pressure and tem- 
perature. The flow was assumed to be l aminar in the first 
case and turbulent in the second. The numerical simulations 
were carried out for approximately 100 (isec, during which 
the projectile would have increased its velocity by approxi- 
mately 20 m/s, assuming a typical acceleration of 20,000 g. 
Since this velocity increase represents less than 1% of the 
projectile’s speed, its acceleration was ignored in the present 
calculations. Also, Bruckner et al [2] demonstrated that the 
acceleration terms in the governing equations can be 
neglected for accelerations less than approximately 20,000 g. 

Case 1 

This case considered the ram accelerator configuration 
shown in Fig. 5a. The explosive gas mixture in the ram 
accelerator section was H 2 + 3 50 2 + 1.5Ar . The projectile’s 

speed was 2136.5 m/s, which corresponded to a Mach num- 
ber of 6.5 in the buffer section and 6.065 in the ram accelera- 
tor section. The flow was assumed to be l aminar , and a 
constant wall temperature of 300 K was specified at the pro- 
jectile surface. A two-block 210x110, 80x159 nonuniform 
grid was utilized. 

The time evolution of the flowfield is shown in Fig. 6 in 
the form of nondimensional temperature 77T contours. 

Fig. 6a shows the projectile just before entering the ram 
accelerator section. The reflected shock wave from the tube 
creates a small separation of the boundary layer, which 
grows slowly with time. At t = 44.58 p.sec (Fig. 6b) igni- 
tion occurs in the projectile boundary layer. Combustion 
spreads both downstream and towards the ram accelerator 
tube (Fig. 6c). A shock-induced combustion wave is estab- 
lished and then reflected from the tube wall, as shown in Figs 
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6d and 6e. The shock-induced combustion wave produces a 
large pressure over the back of the projectile, and a positive 
thrust begins to be generated at t = 69.4 (isec (Fig. 7). At 
t = 68.65 (isec the reflected conical shock created by the 
small ramp ignites the mixture, creating a new shock- 
induced combustion wave (Fig. 6f). This wave magnifies the 
separation in the projectile’s boundary layer, and combustion 
spreads upstream through the boundary layer as shown in 
Figs. 6 g-i. This simulation required 3300 iterations and 11.2 
hrs. of CPU time on a Cray C90 computer, with a maximum 
CFL number of between 3 and 11. 

The net thrust on the projectile is plotted in Fig. 7. Dur- 
ing the transition into the ram accelerator section the drag 
decreases, because the Mach number of the flow ahead of the 
projectile abruptly drops from 6.5 to 6.065. Therefore, as the 
projectile penetrates into the ram accelerator section the 
wave drag decreases progressively. Positive thrust is pro- 
duced after the first shock-induced combustion wave hits the 
projectile. A new peak is formed after the second shock- 
induced combustion wave is established. 

Case 2 

This case considered the ram accelerator configuration 
shown in Fig. 5b. The 30° ramp was shortened, and the pro- 
jectile’s shape modified with the aim of improving perfor- 
mance. The explosive gas mixture in the ram accelerator 
section was + 3.75 + 0.25Ar . The projectile’s veloc- 
ity was 2136.5 m/s, corresponding to a Mach number of 6.5 
in the buffer section and 5.863 in the ram accelerator section. 
The flow was assumed to be turbulent, and a constant wall 
temperature of 300 K was specified at the projectile surface. 
The Baldwin-Lomax turbulence model [19] was used for this 
calculation. A 315x110 nonuniform grid was utilized. 

The time evolution of the fiowfield is shown in Fig.8 in 
the form of nondimensional temperature T/T m contours. 

Fig. 8a shows the projectile at a moment just before bursting 
the diaphragm. The initial flow development is similar to 
that described for the previous case. That is, ignition begins 
in the projectile’s boundary layer (Fig 8b), and a shock- 
induced combustion wave is established (Fig. 8c). This wave 
is then reflected from the ram accelerator tube (Fig. 8d), and 
when it reaches the projectile surface (Fig 8e) a large pres- 
sure is established over the back of the projectile, producing 
positive thrust at t = 58.9 (isec (Fig. 9). In this case the 
reflected shock created by the ramp is not strong enough to 
ignite the mixture. In Fig. 8g combustion is seen to propa- 
gate upstream through the boundary layer. The combustion 
occurring inside the small separation bubble forces it to 
expand. The boundary layer combustion continues to propa- 
gate upstream (Fig. 8h), ultimately resulting in unstart of the 
r am accelerator (Fig. 8i). This simulation required 3400 iter- 
ations and 7.8 hrs. of CPU time on a Cray C90 computer. 


with a maximum CFL number of between 4 and 10. 

The net thrust on the projectile is plotted in Fig. 9. The 
initial development of the thrust force is s imil ar to that for 
the previous case. That is, the drag decreases progressively 
as the projectile enters the ram accelerator section, where the 
speed of sound is higher (and the Mach number therefore 
lower). There is then a sudden jump to positive thrust when 
the reflected shock-induced combustion wave reaches the 
projectile surface. The thrust shows a short increase prior to 
unstart. 

Figures 10 and 11 show the pressure distribution on the 
projectile at various times for the two cases previously 
described. The results for Case 1 are separated into surface 
pressure distribution (Fig. 10a) and projectile base pressure 
distribution (Fig. 10b). At t = 49.92 (isec the surface pres- 
sure is s imil ar to that observed just before diaphragm burst- 
ing (r = 3.27 nsec), except that the pressure level at the 
nose is smaller due to the transition to a lower Mach number 
flow. At t = 76.64 (isec and t = 91.78 (isec the high 
pressure established over the back of the projectile, due to 
the shock-induced combustion wave, can be clearly seen. 
The pressure plots at these times show four peaks. The two 
small peaks are both caused by boundary layer separation. 
The first of the larger two peaks is produced by the ramp and 
the second by the shock-induced combustion wave. The 
pressure distribution on the projectile base is presented in 
Fig. 10b. The base pressure increases significantly during the 
transient phase and then decreases to a level somewhat 
higher than that prior to combustion. 

The pressure distribution for the second case is given in 
Fig. 11. The plots at t = 45.78 (isec and t = 58.09 (isec 
give the distribution after ignition, but before positive thrust 
is being produced. At t = 70.51 (isec positive thrust is 
being created, and the high pressure over the back of the pro- 
jectile is evident in the figure. After unstart 
(r = 89.24 (isec) a large pressure is established over the 
ramp, resulting in a s mall , but negative, total thrust at this 
time. 

Conclusions 

A numerical investigation of the temporal evolution of 
the reacting fiowfield established during projectile transition 
from the launch tube into the ram accelerator section was 
presented. A methodology for simulating the reacting flow 
establishment was described, and computations for two con- 
figurations were presented to illustrate the capability of the 
numerical approach. The efficiency of our time-accurate, 
fully implicit method was demonstrated by computing high- 
speed, reacting, turbulent flows at CFL numbers as high as 
10 . 
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In this study, ignition was always obtained in the bound- 
ary layer, even when the projectile surface was cooled and its 
temperature maintained at 300 K. The combustion in the 
boundary layer spread downstream and towards the ram 
accelerator tube, establishing a shock-induced combustion 
wave. Subsequently, combustion also propagated upstream 
through the separated boundary layer. In one of the cases 
studied, this process unstarted the ram accelerator. 

This woik suggests the need to exert some control of the 
boundary layer flow, for example, through the use of bound- 
ary layer bleed or, as proposed in Ref. [12], a pure hydrogen 
core maintained inside a thin-walled balloon. The latter 
method would not only prevent boundary layer combustion, 
but would reduce the gasdynamic heating of the projectile. 
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Figure 1. Schematic of the University of Washington’s ram accelerator facility [1-3]. 
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Figure 2. Ram accelerator configuration considered in the present study. 
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(a) t = 3.27 nsec 



(b) t = 44.58 



(c) t = 49.92 



(g)t= 76.64 



(h)t= 84.16 



(i) t = 91.78 
Figure 6. continued. 


(d)t= 55.90 



(e) t = 60.94 
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(f) t = 68.65 

Figure 6. Nondimensional temperature T/T contours 
shown reacting flow establishment in ram accelerator. 


Figure 7. Nondimensional net thrust force on ram accel- 
erator projectile. (Here, is the maximum cross-sectional 
area of the projectile.) 


9 











Page intentionally left blank 



11 II M 7J 9 J 11J 



(a) t = 2.53 p.sec 



(b) t = 45.78 



(c)t = 48.33 



(d)t = 53.17 



(e)t= 58.09 




(g)t = 70.51 



(h) t = 77.58 



(i) 1=89.24 
Figure 8. continued. 
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(f)t = 63.04 

Figure 8. Nondimensional temperature 7VT contours 
showing reacting flow establishment in ram accelerator. 


Figure 9. Nondimensional net thrust force on ram accel- 
erator projectile. (Here, is the maximum cross-sectional 
area of the projectile.) 
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Figure 10. Nondimensional pressure distribution at various times; (a) projectile surface; (b) projectile 
Note: vertical distance y is measured from symmetry axis and L (= 10 cm) is a length scale. (Case 1) 
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Figure 11. Nondimensional pressure distribution along projectile surface at various times. L (= 10 cm) is a length scale 
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14. SUBJECT TERMS 

Detonation waves; Reacting flows; TVD schemes 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 


15. NUMBER OF PAGES 

13 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 





















